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ABSTRACT 


The  general  purpose  o^‘  the  work  reported  here  is  to  obtain  the  best 
possible  feasible  signal  processing  algorithms  for  estimating  from  radar 
data  characteristics  and  trajectory  parameters  of  bodies  moving  in  the 
air. 


This  report  starts  with  a--  extensive  set  of  numerical  results  showing 
the  effectiveness  of  the  simplified  dynamics  equations  evolved  Ln  approxi¬ 
mating  pertinent  characteristics  of  trajectories  calculated  by  fche  BBL 
point-mass  model.  Next,  a  mathematical  derivation  of  the  optimal  recursive 
algorithm  for  filtering  and  smoothing  radar  data,  formulaticxi  for  which  was 
presented  in  our  previous  Qiiarterly  Progress  Report  [2],  is  presented.  A 
technique  for  evaluation  of  recursive  filtering-smot ihing  algorithms  by 
deterministic  calculations  is  then  presented.  This  technique  will  enable 
the  assessment  of  expected  random  and  bias  errors  in  trajectory  estimation 
efficiently  so  as  to  avoid  the  need  for  a  large  nuinber  of  simulation  runs. 
In  closing,  plans  are  described  for  presenting  numerical  resixlts  obtained 
in  the  final  report.  These  will  include  true  minimvun  variances  in  trajec¬ 
tory  parameter  estimation  (assiuuing  idealized  processing),  demonstration  by 
simulations  and  deterministic  calculations  of  the  effectiveness  of  the  re¬ 
cursive  smoothing  and  filtering  algorithms  obtained. 
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FOREWORD 

This  report  describes  work  done  frcaa  1  November  1971  through  3I  January 
1972  at  the  Moore  School  of  Electrical  Engineering,  University  of  Pennsylvania, 
under  contract  number  DAAB07“7i“C-oei2  with  U.  S.  Army  Kectronics  Command  for 
research  entitled  "The  Estimation  of  Trajectories  from  Radar  Data".  The  cogni¬ 
zant  technical  personnel  at  USA  ECOM  are  Dr.  Leonard  Hatkin,  head  of  the  Radar 
Technical  Area  and  Mr.  Donald  Foiani,  CSTA  Laboratory,  Evans  Area,  Fort  Monmouth, 
N.J.  07703. 
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TABLES  (continued) 


TABLE  2c  Backtracking  Errors  Due  to  Simplified  Dynamics  Approximations 
for  l^^tm  Shells  . 


TABUS  2d  Backtracki.*  /^  ’Srrors  Due  to  Simplified  Dyn®®ic6  Approximations 
for  0-inch  .  -^fills  . . 


TABIE  3  Ratio  of  Trajectory  Error  to  X-coordinate  Error  A  as  a 


1. 


INTRODUCTION 


As  indicated  in  the  previous  Quai-terly  Progress  Reports  this  work 

is  concerned  with  the  general  signal  processing  problem  of  estimating  charac¬ 
teristics  and  trajectory  parameters  of  a  body  moving  in  the  air  from  noisy 
radar  observations. 

The  main  thrust  during  thir.  quarter  has  been  to  evaluate  the  simplified 
dynamics  equations  evolved,  and  to  complete,  refine,  debug  and  conduct  sample 
runs  of  computer  programs  to  do  the  following  tasks:  a)  obtain  algorithms  to 
perform  optimal  recursive  smoothing  and  filtering  from  simulated  noisy  radar 
dataj  b)  deterministic  evaluation  of  expected  filter-smoother  algorithm  per¬ 
formance}  and  c)  determination  of  idealized  or  minimum  variance  sensitivity 
of  trajectory  parameter  estimation  accuracy  to  radar  observation  errors. 

The  evaluation  of  the  £i’'tplified  Dynamics  Equations  is  presented  in  Sec¬ 
tion  2  to  follow.  Revised  values  for  universal  drag  and  drift  curves  are  pre¬ 
sented,  followed  by  a  set  of  runs  of  a  wide  variety  of  projectile  trajectories 
showing  resultant  backtracking  errors  due  to  the  simplified  dynamics. 

The  detailed  formulation  required  to  perform  optimal  recursive  smoothing 
and  filtering  of  noisy  radar  data  was  presented  in  Section  4  of  ref.  2.  The 
required  mathematical  derivation  of  the  algorithm  is  included  here  in  Section 
3  below.  Details  of  the  simulation  program,  performance  evaluation  for  track¬ 
ing  and  parameter  estimation  of  trajectories  will  be  included  in  the  final  re¬ 
port. 

In  Section  4  is  presented  the  derivation  of  a  deterministic  method  for 
evaluation  of  expected  filter-smoother  algorithm  performance.  This  technique 
can  account  for  bias  errors  as  well  as  random  errors  in  radar  observations, 
for  both  the  optimal  filter-smoother  algorithm  described  in  Section  3  ®nd  also 
arbitrary  non-optimal  recursive  algorithms  that  may  be  employed. 

The  program  for  the  idealized  sensitivity  study  has  been  debugged  and 
results  are  being  prepared  for  the  final  report.  These  results  will  be  use¬ 
ful  in  assessing  the  efficiency  of  the  recursive  algorithms  developed  as  well 
as  providing  tradeoffs  as  to  significance  of  different  types  of  errors. 


The  development  of  simplified  dynamics  equations,  incorporating  approxima¬ 
tions  to  drag  and  drift  accelerations,  was  presented  in  ref.  2,  Section  2.  The 
work  presented  here  consists  of  a  refined  set  of  numerical  constants  needed  for 
the  appro>iraations  used  plus  an  extensive  set  of  evaluation  results  for  the  four 
different  projectiles  considered. 

The  refined  set  of  approximation  constants  is  presented  in  Table  1.  Three 
sets  of  polynoiaial  coefficients  are  given  defining  "universal"  curves  for  drag, 
drift,  and  spin. 


Evaluations  of  the  Simplified  Dynamics  approximations  are  oresented  in 
Tables  2a,  b,  c,  and  d  for  the  105mm,  155®ai#  17^am,  and  8-inch  Projectiles, 
respectively.  In  each  case  resuits  are  shown  for  a  variety  of  charges  and 
quadrant  elevation  angles.  Computations  of  backtracking  accuracy  are  per¬ 
formed  for  reversal  times  given  in  multiples  of  ten  seconds. 

The  principal  factors  affecting  the  errors  are  (i)  time  of  reversal,  corr¬ 
esponding  in  the  use  of  the  system  to  the  time  at  which  the  state  vector  esti¬ 
mate  is  made;  (ii)  quadrant  elevation,  and  (iii)  projectile  type. 

Item  (ii)  is  a  geometric  sensitivity  factor.  For  small  quidrant  elevation 
angle  QE,  the  backtracking  trajectory  is  nearly  parallel  to  the  ground  as  it 
nears  the  launch  point.  A  small  vertical  error  is  magnified  by  the  cosecant  of 
the  angle  Qg.  The  values  of  the  error  in  X#  in  Table  2,  multiplied  by  sin(QE) 
give  values  that  cluster  together  much  better  than  the  values  oC  the  error  as 
given  in  the  Table.  This  is  illustrated  in  Table  3>  showing  the  ratio  of 

the  normal  distance  of  trajectory  computed  by  the  simplified  dyiamics  from  the 
true  trajectory  at  the  launch  point,  to  Z  ,  the  longitudinal  error  in  launch 

A 

point  due  to  the  simplified  dynamics  approximations. 

At  10  second  reversal  time,  the  largest  value  of  (error  in  X)(sin  QE)  is 
4.7  meters,  for  the  155mm  shell  at  QE  =  2l6  mils.  At  20  second  reversal  time 
the  largest  value  is  2^.G  meters,  for  the  same  trajectory.  If  the  155n>n  shells 
are  excluded,  the  largest  error  products  arc  2*5  and  15.6  meters  respectively 
at  JO  and  20  second  reversal  time,  both  for  105njn  shells. 

The  sensitivity  to  projectile  type  suggests  that  the  aerod.mamic  drag  func¬ 
tion  fit  could  be  tailored  better  to  apportion  the  errors  over  the  different 
projectiles,  so  that  the  155mm  ®nd  perhaps  the  105mm  errors  would  be  decreased 
at  the  expense  of  the  175mm  and  8-inch  errors.  For  the  latter  two  projectiles, 
the  largest  values  of  (error  in  x)(sin  QE)  were  2.1  meters  at  10  sec.  reversal 
time  and  12.5  meters  at  20  sec.  reversal  time,  both  for  the  175mm  shell  at 
QE  =  244. 

For  all  the  cases  studied,  the  importance  of  backtracking  zYcxn  a  smoothed 
state  vector  as  soon  as  possible  is  apparent.  The  drag  approxiioaticm  leads  to 
unacceptably  large  errors  for  reversal  times  much  in  excess  of  ten  secemds. 

This  points  out  the  importance  of  the  "smoothing"  as  opposed  to  "filtering" 
features  of  the  recursive  algorithm  being  developed. 
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St^ersedas  Table  B*2  in  Secoad  Qioarterly-  Report 


Table  2b  Backtracking  Errors  Due  to  Simplified  Dynamics  Approximations  for  1^^  ran  Shells 


Table  2d  Backtracking  Errors  Due  to  Simplified  Dynamics  Approximations  for  8-inch  Shells 
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3.0  DEVELOPMENT  OP  ESTIMATION  ALGORITHMS 

3*1  Introduction.  This  section  will  be  devoted  to  the  developoient  of  the  non¬ 
linear  filter-SDoothing  algorithm  presented  in  the  second  quarterly  report.* The 
derivaticm  is  based  on  a  gaussian  approximation  to  optimal  conditional  mean  es¬ 
timation  and  follows  similar  developments  described  by  Wishner  et.al.  [7]  for 
extended  Kalman  filtering.  The  filtering  equations  are  a  modification  of  their 
results;  the  smoothing  equations  are  new.  An  Euler-type  integration  scheme  is 
assumed  for  the  system  dynamics. 


3.2  Problem  statement.  Consider  the  problem  of  determining  sequential  estimates 
of  the  n-dimensial  state  vector  x^  =  of  a  dynamic  system  given  corrupted  dis¬ 
crete  m-dimensional  observations  (or  meisurements ) 

^  =  h(x^)  .  (1) 

!Hie  h(*)  is  an  m- vector  nonlinear  function  of  the  state,  and  the  set  {vj^,te=i,. . .} 
is  a  white  gaussian  sequence  with  r  x  r  covariance  matrix 


I  \ 


V  ''i  *  =  Vkj 


R  I  The  evolution  of  the  state  vector  is  described  by  the  vector  (ito)  stochastic 
p  I  differential  equation 


dx^  =  f(x^,t)dt  +  g(x^,t)dp^ 


where  f(*,t)  is  an  n-vector,  g(*,t)  is  n  X  r,  and  p  is  an  r-vector  Brownian 
motion  process  with 


E{dp^,dp5  =  Qt 


The  i^tial  state  is  assumed  to  be  a  vector  of  gaussian  random  variables  with 

mean  x^  and  covaricuce  matrix  cov(x  ,x  \z  ). 
o  0  0'  o 


3-3  Gaussian  approximation  for  ccmditional  mean  estimaMon.  Ihe  criterion  for 
optimal  estimation  will  be  the  minimization  of  the  conditional  expectation  of 
the  quadratic  form. 

*  I 

Ref.  2,  Section  4.0. 
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XjJj^  Sj  where  S^  is  e«i  arbitrary  positive  definite  weigjhting  matrix.  The 
vector  Xjjj^  =  Xj  ”  *jlk  error  made  in  estimating  the  state  x^  with 
Xjjj^  an  estimate  bused  on  the  given  realization  of  observations 


TvT 


(5) 


The  estimate  that  minimizes  this  criterion  is  known  [U]  to  be  equal  to  the 
expectation  of  x^  conditioned  on  Zj^ 


=  EtxJzJ  . 


*jlk  j'“k^ 


(6) 


lliis  estimate  is  commonly  described  as  the  minimum  variance  or  minimum  mean- 
square  error  estimate. 

Unfortunately  e{x. |z. ;  cannot  be  determined  exactly  for  the  general  non- 

J  K 

linear  problem.  Ohis  is  a  consequence  of  the  fact  that  an  infinite  number  of 
terms  is  required  to  completely  characterize  the  conditicmal  probability 
density  function  p(xj|Zj^)  t43.  If  however  it  can  be  shewn  that  the  and  Zj^ 

are  vectors  of  gaussian  random  variables,  the  density  function  can  be  fully 
described  by  the  following  expressions  for  the  conditional  mean  and  auto  co- 
variance  matrix: 


Etxjlz^)  =  e(xj|z^J  - 


(7) 


cox(Xj.Xj|z^)  =  cov(Xj.Xjla^.^)  -  cov(x^,z^|z^_^)cov-^(i^,z^|z^.^)cov(.^,Xj|z^.j^ 


where 


cov(x^Xplz^)  =  e{{x^  -  E(xjE^.3)(Xp  -  eCx^Iz^DIz^} 


(8) 


(9) 


We  are  motivated  by  the  relative  sin^licity  of  equations  7  6  to  assvuae 

that  for  x  gaussian,  x .  will  be  at  least  nearly  gaussian.  Sorens^  and 
®  0 
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Stubberud  [.4]  explored  this  idea  rather  thoroughly  by  assumirifj  the  a  posteriori 
density  function  to  be  expanded  in  an  Edgeworth  series 


p(x_,|z^)  =  ?(x^K)[l  CjHjCx^lz^)  +  ^  C4fl4(Xj|z^)+  .  .] 


(10) 


wheie  'p{‘ )  is  the  conditional  gaussian  density  function  and  the  H^(* )  eu'e  hermite 

polynomials.  Iheir  experience  with  filtering  algorithms  would  ndicate  that  the 
inc3.usion  of  the  higher  order  terms  in  equation  1.0  is  often  onl.'  of  marginal  value 
and  is  computationally  tsgjroctical  for  all  but  scalar  and  simple  second-order  systems. 
Other  authors  [e.g.  3,5]  have  investigated  dil'ferent  expansions  for  Mith 

similar  results.  Since  our  desire  is  to  develop  a  filter-smoothing  algorithm 
that  may  be  realistically  implemented  for  a  system  relatively  high  order,  we  are 
by  practicality  constrained  to  the  assumption  that  equation  10  nay  be  truncated 
after  the  first  term,  i.e.  that  the  x^  and  are  jointly  gauss.. an.  ihe  experi¬ 
ence  of  other  authors  [summarized  in  reference  8 ]  with  the  reenixy  vehicle  track¬ 
ing  problem  would  indicate  that  this  is  a  reasonable  assumption.  Thus  recalling 
equation  6  we  have  for  the  evolution  of  the  optimal  estimate  due  to  the  observa¬ 
tion  of  2^. 


and 


cov{x  ,xX)  =  covr(x  ,x  Iz^,^)  -  co»(x 


J  3 


where  now 


(32) 


cov(x.,Xj^!z._^) 


(13) 


3*4  Taylor  series  approximations  for  estimate  updating.  Approximations  to  the 
elen.ent  terms  in  aquation  li  and  12  ore  c'otair.ed  by  expnuding  the  various  non¬ 
linear  vector  functions  in  a  Taylor  series  about  .some  nominal  fc-ajectory  state. 

A  logical  choice  for  this  nominal  state  would  be  the  current  best  estimates  of 
Xj^.  Expanding  h(» )  to  first  order  terms  in^^jj^  1  obtain  foi‘  the  measurement 


\  =  *’<\|k-i>  *\*kik-i  -"'k 


(14) 
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a'teM. 


and  thus  taking  expectations 


*k|k-l  ""  ^^\|k-l^  * 

lij^  is  defined  as  the  m  x  n  Jacobian  matrix  of  h(- )  with  element;: 


S  h  (x) 

Wa.^  = 


6  1*  =  '-klk-l 


The  approximate  measurement  residual  z^j-.  ,  nsay  now  be  determined  by  sub- 
trat:ting  equation  l6  from  l4, 


\lk-l  \  \|k-l  ■*■  \ 


Subr.tituting  the  above  ejqp-ession  into  the  defining  term  for  the  measurement 
conditional  covariance 


we  deduce  after  some  raanipiJ.atior.  that 


=  “k“''tV-kl\-i>^  *  \ 


The  following  first-order  approximation  to  the  cross-covariance  matrix 
cov(x^,2j^|z^  1^  determined  similarly  from  equi^tion  17 


C0vfc,2^1r^_l)  =  COv{x 


If  in  equation  20  j  equals  k,  i.e^  only  the  filtering  problem  is  being  con- 
sidcored,  then  equations  19  and  20  would  be  sufficient  and  the  development  of  the 
estimate  update  equations  would  be  complete.  !niat  is,  substituting  equations  19, 
19,  and  20  for  terms  in  U  and  12  one  obtains  for  the  filter 


,.  „  ..  -,- ►o.'-i"*  -•*•  .,-■»-.,«««■■.'«>>  v■■r'^■ w'Vu»*!™rww‘»wn''^»w.mw^rii^«niSKV|nj^^ 


(21) 


\|k  =  ^|k-l  * 


where  we  denote  the  innovation 


\  =  \  -  “(^ik-i* 


witti 


Sk  =  (24) 

It  may  be  noted  the  although  equivalent  to  Wishaer's  [7]  equations  these  ex¬ 
pressions  are  formulated  somewhat  differently.  The  reason  will  become  clear 
irtien  we  determine  the  equations  for  smoothing. 

In  order  to  inclement  fixed-point  smoothing  an  expression  equivalent  to 
equation  19  is  needed  for  updating  the  cross-covariance  matrix  cov(x.,x.  |z.  ). 

Such  on  expression  may  be  deduced  from  the  following  lemma:  ^  ^ 

leaaa  I.  if  u^  and  u^  are  vectors  of  jointly  gaussian  random  variables,  then 

cov(u^,U2lu^)  =  cov(u^,U2)  -  ccv(u^,u^)cov‘^(u^,u^)cov(u^,U2)  (25) 

A  proof  of  Lemma  l  is  included  in  Section  3  below. 

3y  identifying  x^,  and  with  and  u^  respectively  in  equation 

25  and  recalling  that  Zj^  = 

C0»(kj.XklZk>  '  “''('=J'*kl\-l>  -  “''(*J'^i2k.l>“''‘^(V^I\-l>'""(V*klVl> 
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Taken  together  equations  19,  20,  and  26  describe  the  necessary  relationships  for 
updating  the  cross-covariance  terms  cov(Xj,Xj^|Zjj)  and  cov(Xj^,Xj jZj^)(=cov(Xj,Xj^l2^)  ) 

Making  the  appropriate  substitutions  we  determine 


co.(Xj,x^|z^)  =  cov(x^,x^|z^.j)Ll  -  S^H^cov(x^,x^lz^_^)] 


We  shall  use  equation  27  along  with  the  expression 


(27) 


*jlk  =  *jlk..l 


(28) 


•  to  define  the  nonlinear  fixed-point  smoothing  estimate  updating.  !l9ie  smoothing 

I  is  initiated  by  the  filtered  estimates  x^j^  and  cov(Xj,x^(Zj),  or  if  J  =  0,  then 

f  obviously  by  x  and  cov(x  ,x  ).  It  should  be  noted  that  with  the  above  formula- 

\  000 

I  tion,  fixed-poin^  smoothing  introduces  no  additional  matrix  inversion  operations 
I  to  the  estimation  algorithm. 

I, 

Assuming  the  necessary  extrapolation  equations  describing  the  evolution  of 
I  the  estimates  between  samples  are  available,  the  algorithm  is  complete  even 

I  thou^  no  expression  has  been  determined  for  the  auto-covariance  cov(x  ,x  jz^). 

•  We  include  such  an  expression  for  completenessj  J  J  it 


The  development  of  the  dynamic  equations  for  the  various  estimates  and  asso¬ 
ciated  covariance  matrices  between  sample-  will  now  be  considered. 

3.5  Estimate  extrapolation.  Recall  that  the  evolution  of  the  system  state  vector 
is  assumed  to  be  described  by  the  expression 


dx^  =  f (x^,t)dt  +  g(x^,t)d^^ 


(3) 


Proceeding  as  before,  we  expand  f(x^,t)(for  ^  •  first-order  Taylor 

series  about  the  estimate  * 


♦  Pj  +  g(Xj,t)d9^  ^  »k 


(30) 
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E 


where 


a,p 


(31) 


X  =  X 
t» 


tlk-1 


In  li^t  of  our  discussion  on  minimum  variance  estimators  and  equation  6,  we 
determine  the  following  dynamic  equations  for  the  stiboptimal  estimate 


^  *tik-i  -  ‘■(*t|k-i'*>“ 


(32) 


and  estimate  error  between  measvirement  samples 


■=  “‘tlk-l  '  ''t  *tlk.l 


(33) 


In  order  to  facilitate  on-line  estimation,  the  dynamic  equations  will  be 
discretized  by  sinq>le  £uler  numerical  integration.  Integrating  equations  32 
and  33,  we  hav»  for  tj^  *  +  dt 


^|k-l  '  Vljk-l  *  *'^\-llk-l'\-l^ 


(34) 


and 


^jk-1  -  Ll  +  F  k-l^^^Vllk-1  ■*■  ®^*k-l'\-l^Vl 


(35) 


where  the  w^  are  defined  to  be  elements  of  a  white  gaussian  sequence  with 


'’'"j  ’'k’^i '  “  “j-i  “jk  • 


(36) 


In  ordex'  to  sin^lify  the  notation,  we  shall  define  the  following  terms 


Althou^  treating  the  noise  process  in  this  manner  is  certainly  an  abuse  of  the 
theory  of  stochastic  differential  equations  (see  ref  [4j,  chapter  4  for  exanqple}, 
it  will  suffice  in  our  context  of  the  estimation  problem. 
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Errors  introduced  into  the  estimate  extrapolation  by  discretizing  equation  33 
may  be  reduced  somewhat  by  including  second  order  terms  (in  At)  in  the  e]Q)ansiony 
i.e. 


^|k-l  "  Vllk-1  *  ^^Vllk-l»Vl^  *  2  V/^Vllk-r\-l^ 


(40) 


Since  the  inclusion  of  the  additional  factor  does  not  require  that  we  determine 
any  additional  matrix  terms  confuted  anyway],  increased  accuracy 

may  be  obtained  at  little  computational  cost.  As  a  consequence,  equation  Uo 
will  be  esqployed  in  preference  to  equation  34  for  estimate  extrapolation. 

Approximate  equations  for  deter- 

Bineu  by  substituting  equation  39  into  the  appropriate  defining  expression.  Thus 
for 


=  ®^lk-l^|k-lK-l^ 


(41) 


we  h*ive 


*  k-1  ^jk-l  *  ®^*k-l'^-l^k-l^ 

t  *  k-1  \|k-i  *  «<Vi’Vi>*k-i^’^l  W 


m 


X 
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cov(Xj^,Xj^iZj^,^)  =  ^  K-l  Vl>  k-1^ 


•  _T, 


Further  reductions  are  not  possible  without  expanding  8 

tj^  about  the  ^  ^  nominal  trajectory.  Proceeding  with  the  Taylor  series 

expansion^  we  determine  to  first  order  in 

cov(xj^,Xj^|Zj^_^)  =  jj.l  +  «^Vl|k-l'VlK-l 

T. 

[  The  expression  for  the  cross  covariance  term  is  determined  in  a  similar  manner 

E 

i  Viik-i  * 

We  summarize  our  development  by  collecting  the  various  results  into  the 
following  theorem. 

Theorem  I.  The  first-order  filter-smoother  for  the  discretized  nonlinear 
system  described  herein  consists  of  the  following  extrapolation  equations  be¬ 
tween  samples 

I  %ik-i  *  Vi|k-i  *  “  ‘■(vitk-i’Vi*  *  ^  ''k-i  vi|k-i’ Vi>  <'*' 

I  «"(v*kl\-x>  “  ‘k-i”''<Vi'ViK-i>*’^  k-1  ■"  «<Vi|k-i'Vi’Vi 


*^^N-l|k-l'*k-l' 


and  the  following  update  equations  at  the  saasples 


=  \ik-i  *  “''<v*kK-i>\'‘k 
*jlk  =  *slk-l  * 

“’'(»J<*kl\>  =  “’'(’‘j»*ki\.l>tl  -  Sk“k“''<*k’*kl\-X>^  (*■'> 

where 

\ "  *k  -  >''*k|k.i>  (“> 

and 

=k  =  “k’'(«k“''(V\l\.i V  ♦  «k>''' 


^  ^gg^qy*;yy^»p.,„.^  mm„..^^  _ _ _ _ 


i  I 


w.  i 


I 


ts 


m 


I 


3.6  Proof  of  Lemma  i 
Define  the  augmented  vector 


U=^ 


“1 

V 


(47) 


and  sulBtitute  it  into  the  foUoiring  augmented  covariance  matrices 


,  T  T '.  r  cov(u^,u,)  ;  cov(iu,u_) 

cov(u,U)  =  E  1^,  (\  ;»2  )/  =  L  -edi(«^»^)-rcoif«^;,^')j 


(M) 


-1..V  •  ■i(5)-,i .  [  ::::«!j 


(49) 


cov(u2,U)  =  E{u2(Uj^^1  1^^)}  -  [cov(u^,u^);cov(u^,i^)3 


(50) 


These  expressions  may  now  he  substituted  into  the  equation  for  the  conditional 
covariance  cov(UjU|u2)i 


cov(U,u|u^)  =  cov(U,U)  -  cov(U,u^)cov”^(u2,u^)cov(\i2,U) 


(51) 


Performing  the  indicated  matrix  operations  ve  obtain 


CwvCi^.UjlujjTeoVli^.Ugj'ujfJ  "  Lc«f(^,Uj^)' -  «»T(^,T^'Jc5r*^(aj,v^Ic<>»r'>3»'^Jl  re»(»^r>^’-~«»Tri^rOj)«>»”^f«jf«j)eovTuj,\;^)J 

(52) 


The  off-diagonal  term  is  recognized  as  the  desired  result. 


M 


I 


'm 


ii 
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4.0  DETERMINISTIC  COMPUTATIONS  FOR  EVALUATING  FILTER-SMOOTHER  PERFORMANCE 

4. 1  Purpose.  The  conventional  way  to  evaluate  a  RATRAN  signal-processing  filter 
is  by  a  program  of  Monte  Carlo  computation,  as  described  in  the  next  paragraph. 

This  memo  outlines  a  deterministic  method  for  the  same  purpose.  The  deterministic 
method  requires  more  programming.  It  gives  in  return  oetter  insight  into  the  mean¬ 
ing  of  the  results  and  freedom  from  sampling  variances. 

These  results  are  directly  applicable  to  evaluation  of  recursive  filter-smoother 
algorit^  such  as  the  one  developed  in  Section  3  above.  While  expressions  are  in¬ 
cluded  in  the  above  work  for  estimate  covariances,  such  expressions  are  limited  to 
the  specific  conditions  required  for  the  optimalization.  The  formalization  presented 
here  enables  computing  the  effects  of  radar  bias  errors,  and  also  errors  due  to  other 
system  and  environmental  deviations  from  conditions  assumed  in  the  optimalization. 

^•2  C^onventional  method.  In  the  conventional  method,  for  which  the  programming  is 
in  progress,  a  trajectory  is  selected  and  a  number  of  state  vectors 


Xj^  ~  ^^^k^  ^  ~  1, ..•,n 

(53) 

The  noisefree  radar  observation  vectors 

^k  ^^*k^  k  =  l,...,n 

(54) 

are  also  confuted.  These  are  used  repeatedly  with  different  samples  of  radar  noise 
vectors  Vj^,  and  possibly  a  radar  bias  vector  giving 


\ ''k 


(55) 


for  the  radar  observation  vectors.  Uie  z's  are  used  as  inputs  to  the  filter  equa¬ 
tions,  which  conq)ute  the  estimated  state  vectors  jL  and  the  errors  in  the  estimated 
state  vectors  * 


(56) 


For  a  number  of  samples  of  radar  noise,  the  mean  of  the  observed  errors  in  (56)  is 
taken  to  be  the  bias  fi(6x)  of  the  filter  in  the  estimate  of  x(t.  ).  The  covsrience 

K 

of  the  error  is  taken  to  be  the  mean  of  (6x  -  £(6x)(6x  -  E(6x))^,  for  the  noise 
saay>les  used. 

4.3  Outline  of  deterministic  method.  The  deterministic  method,  outlined  next, 
does  not  use  radar  noise  samples.  Instead  it  uses  the  mean  and  the  covariance 
of  the  radar  noise.  The  deterministic  computations  follow  the  propagations  of 
the  means  and  covariances  through  the  filter,  instead  of  the  actual  errors  as  in 
the  Monte  carlo  computations. 


(57) 


(57) 

be  the  true  transition  function,  where  as  before  Xj^  is  the  state  vector  at  time 

tj^  and  now  w^^  is  a  parameter  vector.  The  components  of  w^^  are  thoae  scalar  para- 

ineters  whose  effects  on  the  trajectory  are  of  interest:  wind  velocity  components, 
air  density,  etc.  (in  the  cOTventional  or  Monte  Carlo  computations,  it  is  conven¬ 
ient  to  separate  the  bias  errors  from  the  random  errors}  in  the  deterministic 
method  it  is  convenient  to  discuss  them,  and  even  to  compute  their  effects,  con¬ 
currently.  ) 

The  filter,  in  its  prediction  step,  confutes,  using  an  estimated  parameter 
vector  w. 


^|k-l 


(58) 


Where  g(*,* )  differs  from  f(*,* )  in  RATBAN  because  of  the  use  of  the  singlifled- 
dynamics  equations  in  the  filter,  and  also  because  of  trxmcatlon  errors  in  the 
filter  conqputations.  Let 


'P(x,w)  .  g(x,w)  -  f(x,w) 


(59) 


by  definition,  for  any  (x,w). 
The  filter  then  cooqputes 


^  -  \|k-x  * 


(6o) 


Where  A^  is  the  Kalman  gain  matrix,  here  assumed  to  be  precomputed  —  i.e.,  to 
be  a  function  of  x^^  rather  than  of 

Let  be  the  matrix  with  element  ^  coLuan  j,  where 

h^{* )  is  the  ij^  element  of  the  radar  coordinate  vector  and  is  the  jth 
element  of  the  state  vector  x(tj^),  ao  that 


“(^Ik-i)  -  \  ®Slk-i 


(61) 


-21- 


-- 


tdiere 


H|k.l  •  ^Ik-l  '  \ 

is  the  error  in  the  prediction  step  of  the  filter.  Then  (55  )>  (6o)#  and  (6l)  give, 
approximately, 

•  W^Slk-l 

The  confutation  (58)  may  be  written 


=  g(Vi«Vi>  -  *  ”( Vi'^k-i' 


®’‘k|k-l  “  ’’k-l  ‘4-1  ■*■  “i-l  ‘"k-l 


where  is  the  matrix  that  has  in  row  i,  column  J  the  partial  derivative  of 

the  ith  component  of  g(x,w)  with  respect  to  the  jth  component  of  x,  evaluated  at 
similarly  defined  except  that  the  partial  derivatives  are  with  re¬ 
spect  to  the  confonents  of  W|  and 

Vi  -  Vi'Vi) 

^^k— 1  *  '^k— 1  ”  ^k— 1* 


Substitution  of  (64)  into  (63)  gives 


+ ’'k)  +  (I  -  wxi’k-i  *  4-1  ‘vi> 


*  (I  -  S-i 
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Equation  (65)  hold  equally  well  for  the  Monte  Carlo  coiqputations  as  for  what 
foUoirs  below,  in  the  deterministic  computations,  equation  (65)  is  not  used,  but 
instead  the  equation  that  foUows  trom  taking  the  expected  value  of  both  sides  is 
used: 


E(6xj^)  =  Aj^b(x^)  +  (I  .  +  4:;^  6wj^_^) 

♦  (I  - 

In  (66 ),  it  has  been  assumed  that  w  is  deterministic  and  that  EVj^  *  0.  It  would 

be  easy,  but  I  believe  not  useful  in  RAOSiAN,  to  let  w  have  random  cMiponents. 
Each  symbol  in  (66 )  except  E(6xj^)  and  E(6xj^_^)  represents  a  quantity  that  is 

known  before  the  error  confutations  begin.  Therefore  (66 )  may  be  used,  given  an 
input  value  for  the  mean  initial  error  vector  E(6x  ),  to  confute  the  mean  error 
in  each  x^^.  o 

If  we  write 

i'k  “  ♦  (I  -  WK-i  *  («T) 

for  the  forcing  term  in  (65)  and 

=  (I  - 

for  the  loop  gain,  so  that  (63)  becomes 

=  ^k  *  \  ®  Vi 


(68) 


(69) 


it  follows  that 


cov(6Xj^,6xj^)  =  cov(y^,yj^)  +  ‘^^(^Vl'^^k-l^^ 

If  all  the  terms  in  (67 )  are  deterministic  except  for  v.  ,  and  if  Ev.  =  0, 

K  & 
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and  if  further  the  v.  ,  are  uacorrelated  for  different  values  of  k,  it  foUovs 
that  ® 


cov(yj^,6Xj^_j^)  =  cov(6xj^_^,yj^)  =  0 


Then  (TO)  becomes 


cov(6xj^,6Xj^)  =  cov(Vj^,Vj^)A^  + 


Equations  (66)  and  (72)  give  the  desired  results. 

Equation  (66)  gives  also  the  effects  of  radar  bias  and  other  bias  errors; 
small  bias  errors  of  course  have  no  effect  on  covariances. 

For  smoothing  instead  of  filtering,  an  equation  similar  to(6o),  but  using 
a  different  matrix  in  place  of  A. ,  leads  throu^  identical  steps  to  equations 
similar  to  (66 )  and  (72 )« 


4.4  Explicit  role  of  input  earrors.  Equations  (66)  and  (72 ),  being  recursive, 
do  not  reveal  to  casual  inspection  hew  each  input  error  aj^ects  the  state  vector 
estimate.  For  additional  insight,  and  under  some  conditions  for  an  economy  ia 


computation  at  the  cost  of  some  programing  effort  and  some  storage  space,  there 
may  be  value  in  using  instead  (73  }  and  (7^ )  beloir. 


Equation  (66)  iiqilies 


s(6x^)  =  *  Vk-A-a  *•••* 

idiere 


^k  =  -  *k»k)<Vi  ^  4-i  ®Vi>  (7**)  I 

I 

Equation  (73)  shows  the  contribution  of  the  initial  error  input  to  the  Kalmaa  | 

filter,  and  the  effects  of  the  contributions  trom  each  radar  locdc.  | 


(81) 


cov(6xj^,6x^)  »  cov(6xj^,6Xj^)(u^) 

Finally,  the  mean  error  in  estimated  launch  point,  and  its  covariance,  will 
be  computed  from 


and 


where 


i  “V^o  0  0  0  0  0  0 

(0  -y^o  1  0  0  0  0  0^ 


(82) 

(83) 

(84) 


adjusts  for  the  time  when  the  estimated  altitude  is  zero,  instead  of  t^  when  the 
estimated  time-of-flie^t  is  zero.  In  (84),  z^  are  the  downrange,  up¬ 
ward,  and  lateral  coiqponents  of  launch  velocity. 


I 


‘j.o  CONCLUSIONS 

Tlie  results  of  Section  2  show  the  performance  of  the  simplified  dynamics 
approximations  to  a  variety  of  trajectories  computed  by  the  BRL  point-mass 
model.  In  general  it  is  seen  that  the  approximations  provide  reasonably  accur¬ 
ate  computation  of  launch  i>arameters  provided  that  the  time  from  vliich  back¬ 
tracking  is  initiated  is  sufficiently  early.  Error  variations  are  d\ie  chiefly 
to  differences  between  the  projectiles  and  to  geometric  dependence  on  quadrant 
elevation  angle. 

Section  3  demonstrates  the  mathematical  validity  of  the  optimal  smoothing¬ 
filtering  algorithm  presented  in  the  2nd  Quarterly  Progress  Report. 

Section  4  presents  a  method  for  deterministic  evaluation  of  recursive 
filter-smoother  algorithms,  accounting  for  bias  and  modelling  errors  as  well 
as  random  nuise-like  observation  errors.  This  procedure  can  avoid  the  need 
for  extensive  simulation  runs  for  filter-smoother  evaluation. 

3he  remaining  work  planned  in  the  next  few  months  o?  this  effort  will  be 
mainly  concerned  with  the  application  of  computer  programs  and  techniques  des¬ 
cribed  in  this  and  previous  Quarterly  Progress  Reports.  Specific  results  will 
include: 

(a)  Minimum  variances  of  trajectory  parameter  estimation  errors  caused  by 
errors  in  radar  observations  (assuming  true  optimum  processing).  Sensitivity 
of  results  to  radar-target  geometry  and  radar  system  parameters  such  as  track 
time,  data  rate,  and  S/n  will  be  indicated  as  well  as  an  assessment  of  the 
utility  of  doppler  measurement  for  range-rate. 

(b)  Simulations  of  optimal  filter-smoother  algorithm  performance  for  a 
variety  of  radar-target  geometries.  Evaluations  of  performance  will  be  pre¬ 
sented.  Comparisons  will  be  made  with  true  optimum  (see  (a)  above)  and  sensi¬ 
tivities  will  be  determined  with  respect  to  bias  and  modelling  errors.  While 
some  Monte-Carlo  simulation  runs  will  be  performed,  most  of  the  evaluations 
will  be  based  on  the  deterministic  evaluation  technique  described  in  Section  4 
above. 

(c)  Complete  descriptions  and  listings  of  the  computer  algorithms  and  pro¬ 
grams  plus  a  guide  to  their  use  and  an  evaluation  of  their  xx>ssible  utility  in 
an  operational  system. 


